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Abstract 

The Casimir energy of a massless scalar field is semiclassically given by contri- 
butions due to classical periodic rays. The required subtractions in the spec- 
tral density are determined explicitly The so defined semiclassical Casimir 
energy coincides with that obtained using zeta function regularization in the 
cases studied. Poles in the analytic continuation of zeta function regulariza- 
tion are related to non-universal subtractions in the spectral density. The 
sign of the Casimir energy of a scalar field on a smooth manifold is estimated 
by the sign of the contribution due to the shortest periodic rays only. De- 
manding continuity of the Casimir energy under small deformations of the 
manifold, the method is extended to integrable systems. The Casimir energy 
of a massless scalar field on a manifold with boundaries includes contribu- 
tions due to periodic rays that lie entirely within the boundaries. These 
contributions in general depend on the boundary conditions. Although the 
Casimir energy due to a massless scalar field may be sensitive to the physical 
dimensions of manifolds with boundary its sign can in favorable cases be 
inferred without explicit calculation of the Casimir energy. 

PACS: 12.20.-m, 12.20.Ds, 31.30.Jv 
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1 Introduction 



Classically, the energy of a field is always positive. That need not be so for a 
Casimir energy S, generally thought to have its origin in the vacuum fluctu- 
ations of the field [Q. The possibility of £ being negative can be understood 
from the fact that it is the difference in the infinite zero-point energies of 
the field for two systems. Most calculations in the past were performed by 
directly evaluating this difference of infinite zero-point energies. There are 
many articles and texts that consider Casimir effects in this manner. See for 
instance [2J El El for an overview. 

The ultra-violet divergence of the zero-point energy in general reflects 
local properties of the system. From a path-integral point of view the di- 
vergence is due to contributions from arbitrarily short paths that begin 
and end at the same point. These zero- length paths probe the local ra- 
dius of curvature or, for paths touching boundaries, local properties of the 
boundary [3 El El IH1 IS] • In a few favorable situations, ultra-short paths do not 
contribute to the difference of zero-point energies. This in particular is the 
case for rigid disjoint boundaries that are moved relative to each other [TU]. 
One sometimes also considers idealized boundaries whose local deformation 
does not cost energy. An example of the latter are smooth, perfectly con- 
ducting metallic surfaces of vanishing thickness in three dimensions. The 
(local) surface tension of such an ideal surface vanishes |12t lT3j and no energy 
is required to deform it locally. 

A measureable Casimir energy should not be extremely sensitive to the 
(sometimes implicit) ultra-violet cutoff and should depend on global char- 
acteristics of the system only. It otherwise is difficult to disentangle the 
energy required to change the system as a whole from purely local effects, 
for instance due to changes in the local curvature of the space or in the 
transmission of a boundary. 

There thus either is no contribution from ultra-short paths to a mea- 
surable Casimir effect or it must be possible to unambiguously isolate these 
local contributions to the vacuum energy. Most calculations of the Casimir 
energy based on spectral properties and Green function methods do not sep- 
arate length scales explicitly. The regularization and subsequent subtraction 
of divergent contributions often are difficult to motivate physically and it 
is not always apparent how Casimir energies of different systems can to be 
compared. 

The semiclassical evaluation of Casimir energies advocated in^3] relies 
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on an ah initio separation of scales. The approach separates the semiclassi- 
cal contribution to the Casimir energy due to quadratic fluctuations about 
classical periodic rays (paths) from all others. Classical periodic rays are of 
(finite) extremal length and give a semiclassical approximation to the Casimir 
energy that depends on global characteristics of the system only. This part of 
the vacuum energy is naturally finite and does not include ultra-violet contri- 
butions from length scales that are much smaller than the shortest classical 
periodic ray. This is one of the principal conceptual differences to the "op- 
tical" approximation to Casimir energies [T5]. The latter (in principle) takes 
all closed classical paths (not just periodic ones) into account 1 . Closed paths 
can be ultra-short in the vicinity of surfaces and lead to divergent Casimir 
self-energies. The optical approximation therefore has mainly been used to 
obtain a numerical estimate of the interaction energy for rigid bodies. 

It was argued in that a semiclassical evaluation of the Casimir energy 
is often particularly simple and gives the leading asymptotic behavior when 
the Casimir energy is large. This is the experimentally most accessible region 
of parameter space [T7j. 

However, the desired separation of length scales may not always be pos- 
sible: changing the radius of a spherical shell invariably changes the local 
curvature as well. The energy required to achieve a change in radius in this 
case will include a possibly divergent contribution from the local change in 
curvature. 

This suggests dividing systems into classes: the difference in vacuum 
energy of any two systems within the same class being finite. It would re- 
quire an infinite amount of energy to compare systems belonging to different 
classes. Within a particular class, the finite Casimir energy has the univer- 
sal interpretation of a vacuum energy: differences in Casimir energy are the 
finite differences in vacuum energy. 

The spectral density p(E; . . .) is assumed to be a well-defined quantity 
for any system (at least for "free", non-interacting fields). The ellipsis here 
stand for the space Ai, the types of field, the boundary conditions that are 
satisfied and any other qualifiers of the system. For systems A and B of the 

1 Although this appears to be an improvement over the semiclassical treatment, the 
optical approximation to Casimir energies also only includes quadratic fluctuations about 
classical rays. The optical approach in principle could provide a more uniform approxima- 
tion in some cases (but not in all) but inherently is no more accurate than the semiclassical 
approach. An objective comparison of the two methods|16| is rather difficult due to the 
numerical limitations and approximations of this approach. 
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same class, the difference of spectral densities, 



p(E;A-B):=p(E;A)-p(E;B) , (1) 
by definition has a finite first moment, 



— oo < £a_ 



A—B 



i r°° 

- / p(E; A - B)EdE = E vac (A) - £ vac (£) < oo . (2) 
* Jo 



&A-B could be called the Casimir energy of system A with respect to system 
B. There evidently are many equivalent definitions of the Casimir energy of 
a system within a particular class - they are distinguished by the spectral 
density used as reference. The Casimir energy determined by two such sub- 
traction schemes, differs only by a finite amount that is the same for any 
system of a class. Such subtraction schemes are equivalent in all physical 
respects. 

The semi classical Casimir energy (SCE) is defined by a particular sub- 
traction po(E) in each class. I will take advantage of the fact that the semi- 
classical spectral density p(E) is the sum of a part p(E) determined by con- 
tributions from periodic rays and a (often classical) remainder |1 8 j po{E), 

p(E; . . .) = p(E; ...)- p Q (E] class) = -- lim Im g(E + ie; . . .) . (3) 

Here g(E) is the part of the response function due to classical periodic rays. 
For a scalar field the remainder po(E) at least includes the Weyl contribution 
to the spectral density proportional to the volume of M.. In addition po(E) 
may depend on the type of field, the curvature, boundaries as well as other 
characteristics El EES] • 

The SCE S c then is defined as, 

PCO T? 1 /■OO T? 

£JM)= -p{E)dE = lim/ -lmg{E + ie)dE . (4) 

J 2 7T e^O J Q 2 

Since the length of a periodic ray is finite, this contribution to the vacuum 
energy is free of ultra-violet divergences and in general is finite[TQ. The SCE 
of Eq. (J3J may be taken to (at least approximately) represent the vacuum 
energy within a class of systems for which the subtracted spectral density 
Po(E) is the same. It may happen (see the example of the Laplace-Beltrami 
operator on a half-sphere of Appendix A) that a particular class has just 
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one member. The subtraction po(E) in this case is not universal to several 
(> 2) systems. The finite Casimir energy one extracts in this case is peculiar 
to a particular system and is physically quite irrelevant: any change in the 
system requires infinite energy. In several examples studied in Appendix A 
such non-universal subtractions are associated with poles in zeta function 
regularization. The SCE of Eq. on the other hand coincided with the 
Casimir energy of zeta function regularization in all systems I studied for 
which the subtraction has a more universal meaning. 

Although Eq. (J3J) does not directly refer to po(E), this implicit subtraction 
in the spectral density determines the class of systems and thus, in effect, the 
usefulness of the SCE. Other approaches, such as zeta function regulariza- 
tion often give finite answers without specifying what has been subtracted. 
Still other approaches, such as heat kernel expansion, subtract terms whose 
physical implications are not entirely clear jHJ E] and the question whether one 
gains or looses vacuum energy by transforming an elongated ellipsoid into a 
sphere is difficult to answer. To escape this conundrum in the interpretation 
of a Casimir energy, Power [TT?] long ago considered a large rectangular box 
with a moveable wall to define the original Casimir energyjTj for two paral- 
lel conducting plates unambiguously. He in effect was considering a class of 
systems that all have the same total volume, total surface area, edge length 
and number of corners. We will see in Section 14.21 that the implicitly sub- 
tracted spectral density po(E) for a three-dimensional parallelepiped in fact 
only depends on these characteristics. p(E) of a parallelepiped can again be 
expressed in terms of periodic orbits |2"U1 |2"T) 122"] only. 

Restricting the validity of a Casimir energy to a certain class of spaces 
for which the same subtraction in the spectral density gives a finite Casimir 
energy in this sense generalizes Power's procedure to slightly less obvious 
situations. As the example in Appendix A of a massless scalar field on 
S4 demonstrates, (universal) subtractions can go beyond Weyl terms and 
for instance include contributions proportional to the integral of the (local) 
curvature over the whole space. 

As emphasized in a perturbative setting by Barton[23], the physical in- 
terpretation of a Casimir energy depends almost entirely on the (implicit) 
subtraction. This is readily illustrated by a spherical cavity in three dimen- 
sions. The significance [21] of the electromagnetic Casimir energy (which was 
found to decrease with the radius of the cavitv[TT]). relies on the fact that 
this Casimir energy actually determines the physical pressure on the spher- 
ical surface of the cavity. This conclusion is possible only if the (implicit) 
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subtractions in po(E) do not depend on the surface area of the boundary. 
The finite Casimir energy otherwise could only be used to obtain the vac- 
uum energy difference between cavities of the same surface area and would 
not determine the pressure on the cavity surface. That a subtraction propor- 
tional to the surface area is not required in the electromagnetic case, is due 
to the ideal metallic boundary conditions 113] . The situation is less fa- 
vorable for a scalar field [8 J satisfying Dirichlet boundary conditions on such 
a spherical surface. The non-universality of the required subtraction was 
emphasized in[7j. 

Defining the Casimir energy in terms of contributions due to periodic 
orbits rather than by any other subtraction of the spectral density has the 
advantage that this finite part of the vacuum energy may often be evaluated 
approximately. This is of practical use in situations where the exact spectrum 
is not, or is only numerically, known. A rather crude approximation will give 
an estimate of the sign of the SCE in Eq. (j3J) without detailed knowledge of 
the periodic rays themselves. 

The sign of Casimir energies is one of its many puzzles. Without explicit 
calculation, determining the sign of the difference of two divergent vacuum 
energies in general is quite hopeless. Obtaining the sign of the SCE on the 
other hand is much more promising due to the geometrical nature of this 
definition. The overall phase of the contribution to the response function 
from a particular periodic ray is given by a topological winding number j2S]- 

1 will argue that the sign of the SCE can often already be inferred from the 
shortest periodic rays that contribute. 

I first illustrate the approach for single valued (bosonic) fields on smooth 
ci-dimensional manifolds without boundary such as Sd and T^. I then general- 
ize to manifolds with boundaries on which the bosonic field satisfies Dirichlet 
or Neumann conditions. Several examples show that classical periodic rays 
within the boundary must also be considered. In general the contribution of 
these rays to the SCE depends on the boundary condition. When the bound- 
ary is not smooth, as for a parallelepiped, contributions due to periodic rays 
in even lower dimensional spaces have to be included as well. 

2 General Spaces without Boundary 

The conceptually simplest Casimir energy probably is that due to a massless 
single valued bosonic field on a smooth d + 1-dimensional Riemannian space- 
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time without boundary. I will assume that the metric is static in a particular 
frame, i.e. that it makes sense to speak of a (^-dimensional spatial manifold 
A4 and of the energy of a particle. Periodic rays follow geodesies on M. that 
close on themselves. The classical action for a periodic ray 7 then is, 

S~ — (p p ■ dx = p(E)L~, (5) 



where r 7 = dS 7 /dE is the time for the ray to return to its starting point 
on M.. For a massless particle moving at the speed of light, r 7 = L 7 /c and 
thus p(E) = E/c. Note that for periodic rays r 7 is an integer multiple of the 
primitive period t 7 . 

The contribution g(E) of isolated periodic rays to the response function 
is of the form[18 , 

g(E) = ^J2 Ayt^i^-^ 2 . (6) 
7 

In Eq. © the amplitude A 1 is determined by the monodromy matrix associ- 
ated with the ray 7. It is a geometric quantity that (for massless particles in 
vacuum) does not depend on their energy E. v4 7 furthermore is positive and 
real by definition. The integer o~ 7 > is the Maslov-like index of the stable 
and unstable manifolds of the periodic ray[22]. Important for us is that this 
index is a topological winding number. As such it is an additive integer that 
scales directly with the number of times an orbit is iterated. For practical 
calculations it will be useful that 07 may be written as the sum, 

(X 7 = yU 7 + V-y , (7) 

of the number of conjugate points /i 7 between the initial point x and the 
final point x' = x of the periodic ray and of an integer z/ 7 associated with 
the stability of the periodic orbit. 

/i 7 gives the total phase retardation /i 7 7r/2 due to conjugate points (for 
manifolds without boundaries) encountered by the periodic ray. z/ 7 in Eq. ((7|) 
is the number of negative eigenvalues of the matrix W of second variations 
of L 7 with respect to a change of the initial (=final) point of the periodic 
ray HE! , 

5L 7 (x + 6y, x' + 5y) | x=x , = 5y T ■ W(x) ■ 8y . (8) 

Since one of the eigenvalues of the dxd matrix W always vanishes, < v 1 < 
d-1. 
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If all periodic rays are isolated, one can insert Eq. (jUJ) in the definition of 
Eq. (JH). Upon performing the energy integral, the SCE due to only isolated 
periodic rays is of the form, 



The sign of the contribution of a particular periodic ray to the SCE is de- 
termined by the integer a 7 . Remarkably, periodic rays with odd ct 7 do not 
contribute to the Casimir energy. Although an expression like Eq. Q is 
valid only for isolated periodic rays, it can also be used to estimate the 
sign of the Casimir energy of integrable systems. The expression of Berry 
and Tabor |2E] for the spectral density of an integrable system in terms of 
periodic rays is more appropriate (see below), but integrable systems are 
singular in the sense that small deformations of the manifold M. destroy the 
symmetries and generically result in isolated periodic rays. The expression 
of Eq. © is robust in the sense that the SCE changes continuously and in 
particular generally does not change sign if the deformation is small enough. 
Slightly deforming M. to isolate the orbits therefore should allow us to obtain 
the sign of £ c from Eq. © even for integrable systems. In support of this 
conjecture note that the integer cx 7 of an individual periodic ray in Eq. © 
is a winding number that changes only when a new caustic appears or the 
stability of the periodic ray changes. If the contribution of a periodic ray to 
the Casimir energy does not vanish, its sign should not change for sufficiently 
small deformations of the manifold. 

A unique determination of the sign of S c is possible when cos(<t 7 7t/2) 
does not depend on the periodic ray 7. In less favorable situations I resort 
to finding the sign of the contribution due to the shortest periodic rays to 
Eq. ©. The contribution in Eq. Q of a periodic ray that winds n times 
about the geodesic generally decreases in magnitude as 1/n 2 and in some 
cases decreases even faster. If the contribution from primitive periodic rays 
dominates the SCE in Eq. 0, I will estimate its sign by that of the shortest 
primitive periodic rays. For some spaces (see for instance S 2n below), short 
primitive rays do not contribute to the Casimir energy at all or may give 
contributions of either sign. The overall sign of the Casimir energy in this 
case is ambiguous and this estimate fails. One nevertheless might expect S c 
to be rather small in magnitude in such situations and I will write S c ~ 
when the sign cannot be determined from the shortest periodic rays. 
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2.1 (i-dimensional Tori and Spheres 

Obtaining the sign of the SCE is straightforward for a massless scalar on a 
(i-dimensional torus Td = Si x Si x . . . x Si. The curvature of Td vanishes and 
it is a space without boundary. The subtracted spectral density therefore is 
the Weyl-term proportional to the volume of Td only. However, due to the 
translational symmetries this is an integrable system and classical periodic 
rays are not isolated. To estimate the sign of the SCE of a torus using 
Eq. ©, one has to deform it slightly. Such a deformation generally destroys 
all symmetries and gives isolated periodic rays, the shortest of which resemble 
periodic rays of the original torus on its shortest cycles. If the curvature 
remains sufficiently small on the deformed torus, the number of conjugate 
points along a primitive orbit continues to vanish. The length of the shortest 
periodic rays furthermore is a minimum by definition. One thus obtains 
= for the shortest periodic rays of a slightly deformed torus. 
They all give a negative contribution to the Casimir energy in Eq. Q. Since 
this sign does not depend on the particular deformation of the torus, one can 
be confident that, 

£ c (T d ) < for all d= 1,2,... (10) 

This sign is in agreement with that of the Casimir energy due to a massless 
scalar field satisfying periodic boundary conditions on the hyper-surface of 
any (i-dimensional parallelepiped obtained by explicit calculation [22\. 

Eq. Q indicates that the SCE may in principle be of either sign for 
manifolds without boundaries. A non-trivial example is the Casimir energy 
of a massless scalar field confined to a spherical shell Sd of dimension d and 
radius R. Periodic rays follow great circles of radius R on Sd- They again are 
not isolated and the d — 1 dimensional cross section of a bundle of initially 
parallel geodesies is reduced to a point at two anti-podes. Any starting 
point on a geodesic of Sd also is a conjugate point of order d — 1, i.e. it 
is self-conjugate [27| . One can avoid the associated complications by slightly 
deforming Sd in a generic fashion. The primitive periodic rays of the deformed 
sphere should still resemble original geodesies on Sd, but the starting point 
generally no longer is self-conjugate. A pencil of rays emanating from a 
point on a geodesic also generally no longer meets in a single focal point. 
For a sufficiently small deformation, focal points of the sphere will have 
been resolved into a series of d — 1 closely spaced conjugate points of first 
order. On the shortest primitive periodic ray, there generically are only d—1 
conjugate points, the next bunch of d—1 conjugate points occurring just after 
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completion of a full revolution. The reason is that the curvature along the 
shortest primitive ray in general is slightly too small to lead to more than one 
intersection 2 - as for a periodic ray about the waist of a slightly elongated 
ellipsoid. The shortest periodic rays furthermore are of minimal length and 
therefore are stable and v 1 = in this case. The total phase retardation 
of the shortest primitive rays thus is <t 7 7t/2 = /i 7 7r/2 = (d — l)vr/2. Their 
contribution to the Casimir energy of a (deformed) sphere in Eq. © vanishes 
in even dimensions and is of alternating sign in odd dimensions. Assuming 
that the contribution from the shortest periodic rays dominate in Eq. ((SJ), 
one thus estimates that, 



The explicit calculations for d < 4 of Appendix A confirm this not very intu- 
itive pattern for the sign of the Casimir energy of a massless scalar field on 
low- dimensional spheres. A little surprisingly, the Casimir energy vanishes 
exactly for S2, S4 and in fact any S even . Since our somewhat crude estima- 
tion only takes contributions from the shortest primitive rays into account, 
it cannot in itself predict a vanishing Casimir energy. When the contribu- 
tion from the shortest periodic rays vanishes, determining the overall sign 
of the Casimir energy becomes much more involved. To conclude that the 
Casimir energy vanishes one has to show that it is positive for some (small) 
deformations of the manifold and negative for others. Although this can be 
shown for Seven, the argument is no longer "simple" and I will not pursue it 
any further. 

It is amusing to consider further stretching the sphere to an elongated 
cigar-like shape and eventually a long cylinder with end caps. The subtrac- 
tions have to remain constant during this deformation. In two and three di- 
mensions, this amounts to keeping the total surface area, respectively volume, 
constant during the deformation. In higher dimensions, certain moments of 
the curvature and other global characteristics also must not change. As the 
curvature decreases, the d — 1 conjugate points on the shortest (and thus 
stable) periodic rays move beyond the end of the primitive orbit when the 

2 The geodesic distance r/(s) between two nearby geodesies satisfies the linear second 
order equation d 2 rj(s)/ds 2 — —n(s)r/(s), where s is the arc length along the ray and k(s) 
is the Gaussian curvature at s. If k(s) < Air 2 /!* 2 two geodesies at most meet once. This 
condition generally holds for the shortest primitive ray of a slightly deformed sphere. 




< 0, for d = 1 mod 4 
~ 0, for even d 
> 0, for d = 3 mod 4 



(11) 
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sphere is stretched sufficiently. a 1 — > for these orbits and they eventually 
dominate and lead to a negative Casimir energy for a (slightly deformed) 
long cylindrical d- dimensional surface (for essentially the same reason as the 
short rays on a torus). Since the length of the shortest periodic rays on a very 
elongated cylinder is much less than on a sphere of the same ^-dimensional 
volume, the Casimir energy furthermore increases in magnitude during this 
deformation 3 . A free massless scalar field thus would collapse to a fila- 
ment. 

3 Integrable Systems 

A quantitative comparison with the Casimir energy obtained by other meth- 
ods generally is possible only for integrable systems, for which the true spec- 
trum is explicitly or implicitly }28| known. 

In the examples studied in Appendix A the SCE coincides with that 
obtained by zeta function regularization when the analytic continuation is 
uniquely defined. When poles arise and zeta function regularization gives 
ambiguous answers, the associated subtraction in the spectral density is not 
universal either. Universal subtractions on the other hand may include terms 
in po(E) that depend on global characteristics of the curvature [as for the 
Casimir energy on S4 in Appendix A]. 

All our estimates of the sign of the SCE rely on the dispersion relation 
p(E) = E/c of a massless particle, since the integral in Eq. (j3J) had to be 
performed explicitly to arrive at Eq. The sign of the SCE in general 
will differ for other dispersion relations, as for instance for the spectrum of 
the Laplace-Beltrami operator without curvature correction. As argued in 
Appendix A the latter is not the spectrum of a massless scalar HE] . Using 
zeta function regularization, the Casimir energy in this case is negative for 
all Sd,d < 4 30J, but the Casimir energy of S3 is ambiguous due to a pole 
contribution. 

The explicit calculations of Appendix A strongly suggest that the SCE of 
a massless scalar on vanishes exactly for even dimension d: the integrands 
in this case are polynomials in (ER) 2 only. The subtraction po(E) for S eyen 
otherwise would not be universal. The integral over the energy of the response 

3 It would be erroneous to compare the Casimir energy densities of a cylindrical surface 
and a spherical one of the same radius when the subtractions are not the same. 
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function (after Wick rotation) has vanishing imaginary part and there is no 
contribution to the Casimir energy from any periodic ray on S eyen . 

Although this supports the previous estimate of the sign of the SCE for 
low-dimensional spheres, Eq. Q requires periodic rays that are isolated 
and thus cannot be directly applied to integrable systems. Let us therefore 
obtain an expression for the SCE that can be used to determine its sign in 
the (rather special) case of an integrable system. 

Act ion- angle variables are the canonical phase- space coordinates of inte- 
grable systems. The SCE for a d- dimensional integrable system is found by 
applying Poisson's formula (HI 



Here H(T) is the classical hamiltonian expressed in terms of the actions I = 
(Ji, lit ■ ■ ■ , Id) and m is a <i-dimensional vector of integers. The summation 
in Eq. (fT2|) extends over all such vectors except m = = (0, 0, . . . , 0). This 
contribution has been subtracted by the "classical" spectral density, 



Eq. (|T2*j) expresses the SCE of an integrable system as a sum over classical 
periodic trajectories on the invariant torii. The classical action of a trajectory 
with winding numbers m about each of the cycles of the invariant torus is 
S(m) = 27rm • I. The correction to the classical action proportional to 
h is linear in the winding numbers m. It is a topological quantity that 
determines how periodic orbits on the invariant torus are projected onto 
physical coordinate space jSI]. The vector j3 of Keller-Maslov[33 indices gives 
the phase loss from caustics on a periodic orbit (see below). Note that (3 is 
a geometrical quantity that does not depend on I. 

The integrals over the actions of Eq. (|T2*|) are evaluated semiclassically 
at stationary points I(E, m) of the classical action on the energy surface - 
the vector m is normal to the energy surface H(T) = E at I = I(E, m). For 
given energy E, the d — 1 integrations along the (compact) energy surface are 
performed semiclassically by choosing a local frame of actions at I for which 
one axis, say that of 1%, is in the direction of m and all others are tangent 
to the energy surface. Care must be taken with zero modes of the matrix of 



/ integrable \ = J_ ( h(T) e 
c V system J 2h d ^ J U 6 



27rim-[I//l-/3/4] d j 



(12) 




(13) 
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second derivatives 



Hij = PH/dlidlj, i, j = 2, . . . , d . (14) 

The integral over the u dimensional subspace of the zero modes of is 
stationary only when the corresponding frequencies vanish. The sum over 
m is thereby reduced to one over the (d — z/ ) -dimensional vectors n that 
are orthogonal to the zero modes. Denoting the uq- dimensional volume of 
this classical moduli-space by V UQ (E, n), the semiclassical result for the d—1 
integrations along the energy surface in Eq. (|T2jl is, 



integrable \ ^ f°° EdE V Uo {E, n) e ^("--^+)/4 



/ integrable \ f 
V system J £f j 



e 27rin.[I/ft-/3/4] 



,,..„— 2ujhd y/\n/hLu\ d - l -»°\det' H i:j \ 

(15) 

The primed determinant here means that the uq zero-modes have been omit- 
ted in its calculation. v + and v- are the number of positive and nega- 
tive eigenvalues of Hy, u + + i/_ + u = d — 1. The frequency uj = \u\ = 
=i(E,m) f° r a rnassless particle in fact does not depend on its energy. 
The energy dependence of the integrand in Eq. (fTK|) is made explicit by 
noting that the classical action of a rnassless field is S(E, n) = |27rn ■ I| = 
EL n /c where L n > is the length of the periodic ray. For dimensional 
reasons, det' (E) scales as (E / oj 2 ) l ~ d+v ° and V U0 (E) scales as (E/u) uo . 
The integration over the energy E in Eq. (JT5j) then gives the SCE of an 
integrable system without boundaries in the form, 

/integrable\ = f c X* 1 ^ cos{<(5 _ Vq _ v _ )j2) . (16) 

V system J ^ \L n J 

To perform the integral in Eq. (|15j) . the integer 

5- n = (/3 - vo - V-)ji , (17) 

must not depend on E. Since geodesies do not depend on the energy of a 
ray, the number of non-positive eigenvalues of Hij does not depend on E for 
rnassless particles. The Maslov-Keller index /3 n = n ■ f3 depends only on n 
by construction. Note that any constant angle variable of a closed geodesic 
implies that the hamiltonian does not depend on the conjugate action. The 
dimension < u < d — 1 of the space of zero modes of can thus often 
be obtained by inspection. 
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The amplitude in Eq. (fTtjj) is positive by construction and I have 
absorbed all dependence on the scale |n| in the length L n of a class of periodic 
rays. thus is a function of the dimensions of the integrable system (such 
as the volume of the space) that do not scale with the length of the periodic 
orbit. The leading contribution to the Casimir energy of an integrable system 
in Eq. generally is from the shortest rays that contribute to Eq. 
The sign of this contribution depends on the Keller-Maslov index (3 n of the 
rays as well as the number of non-positive eigenvalues of Hy. The sign of 
the contribution from a particular class n of periodic rays of an integrable 
system is again given by, 



with cr n defined by Eq. (fT7j) . The SCE of Eq. (fTT)|) for an integrable system 
thus is of a remarkably similar form as the one in Eq. Q for a system with 
only isolated periodic rays. It is tempting to identify a n in Eq. (|T7|) with 
cr 7 of Eq. (J7|). However, while both <r 7 and a n are topological quantities, <r 7 
in Gutzwiller's expression for the contribution to the response function of 
isolated rays is a topological property of the ray [25] > whereas a n in Eq. (fT7j) 
is a property of a whole continuous family n of periodic rays. The two 
integers depend differently on the winding number of the periodic rays. One 
nevertheless can argue that they coincide for the shortest (and thus primitive) 
periodic rays. 

The Keller-Maslov index /3 n of an integrable system is given by the num- 
ber of caustics a periodic ray encounters (for manifolds without boundary). 
These caustics are created by the family of rays it is a member of. f3 n thus is 
a topological property of the family of rays that generally does not equal /i 7 , 
the number of conjugate points on a single ray of that family. For one, f3 n 
does not depend on the starting point of the ray and is a well-defined integer 
even when this point happens to be self-conjugate (as in the case of spheres). 
f3-a is proportional to the number of times the periodic rays wind about the 
closed geodesic. The number of positive eigenvalues of on the other hand 
is a statement about the curvature at a particular point on the energy sur- 
face determined by the direction of its normal n. Neither this point nor the 
energy surface change with the magnitude of n. (z/ + v_) = (d — 1 — u + ) 
thus does not depend on the winding number of the periodic ray. For cr 7 to 
coincide with a n for all periodic rays 7 as the integrable system is slightly 
deformed, the number of non-positive eigenvalues of has to vanish (as for 
a torus). However, the sign of the Casimir energy in general does not change 




(18) 
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under small deformations, if the phase of the shortest periodic rays remains 
the same. This observation gives the desired connection between cx 7 and a n : 
the two have to coincide for the shortest periodic rays of an integrable system 
and its (sufficiently small) deformation. Below I show that this is the case 
for the previous example of spheres and tori. 

3.1 Covering Spaces 

The previous arguments imply that the sign of the Casimir energy can often 
be inferred from the Keller-Maslov index /3 7 of periodic rays of integrable 
systems. [7 here is some representative of the class n.] It thus is important to 
have a reliable and transparent determination of this index. Keller pTT] gives 
a geometrical construction that generalizes to manifolds with boundaries. I 
recall points of the construction that are relevant here. 

A solution S(q,t) of the Hamilton- Jacobi equation H(q,p = VS,t) = 
E for constant energy E generally is multiply valued and so may be the 
momentum p = VS* itself. At any point x the momentum p = VS* of a 
solution only has a finite number of branches, say m of them. One constructs 
an m-sheeted covering space on which p = VS is a single- valued function by 
associating each of the branches of VS with a separate sheet. Any two 
different sheets % and j are joined together on sub-manifolds on which the 
momenta coincide p\ { = p\-. Such a sub-manifold generically is a caustic 

of the rays or a boundary in the original space. [If VS is defined on only 
part of x-space then only this part is covered.] The advantage of Keller's 
construction is that p = VS becomes a single- valued function on the covering 
space. A family n of periodic rays does not intersect on the covering space. 

The semi-classical Casimir energy of the integrable system is obtained by 
considering the periodic rays on this covering space. The phase is retarded by 
rmi/2 whenever a ray crosses a caustic of m th order in passing from one sheet 
to another. The positive integer m is the number of dimensions by which 
the cross-section of a tube of nearby trajectories is reduced at the caustic. 
Keller's construction does not of itself provide the order of a caustic. The 
latter must be inferred from the behavior of a bundle of nearby rays. 

An example is the construction of this covering space for geodesic motion 
on a (^-dimensional sphere Sd- The SO(d + 1) symmetry of Sd implies that 
geodesies lie in a hyperplane of the R^+i it is embedded in. The vector or- 
thogonal to this hyperplane makes a certain angle 9 with the vertical which 
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is the inclination of the geodesic. It can be chosen as one of the angles of 
the action-angle variables and geodesies of Sd thus fall into distinct classes 
characterized by their inclination. A family of periodic rays of the same incli- 
nation covers an annulus of Sd that is bounded by two Sd-i hyper-surfaces. 
These two d — 1-dimensional caustics are around the polar regions of Sd- Ev- 
ery closed geodesic in the family of solutions with given inclination touches 
the "upper" and "lower" caustic once. One constructs the covering space 
by joining two sheets of this annulus of Sd at the two caustics. A periodic 
ray of the given inclination passes from one sheet to the other every time 
it crosses one of these caustics. In this covering space a family of rays of 
fixed inclination does not intersect. The number of times a periodic ray of 
fixed inclination winds about the annulus also is the number of times the 
periodic ray passes through both caustics. These caustics are of order d — 1, 
because the cross section of a d — 1-dimensional tube of geodesies vanishes 
when they intersect a given periodic ray at the caustic. The phase loss of 
a periodic ray that winds about the sphere n times and crosses In caustics 
is 2n(d — l)7r/2 and thus /3 7 = 2n(d — 1) for any family of periodic rays of 
fixed inclination. The hamiltonian depends on the action conjugate to the 
angle describing the motion on a great circle only (the magnitude of angular 
momentum in the 2-dimensional case). Therefore uq = d — 1 and V- = in 
this case leading to a = (2n — l)(d — 1). For even dimensions d, a thus is 
odd for all n and the contribution to the SCE of all periodic rays according 
to Eq. (|18|) vanishes. This agrees with our previous estimation of the signs 
in Eq. (jllj) based on Eq. © and proves that the SCE of even- dimensional 
spheres indeed vanishes. Appendix A presents explicit results for < d < 4. 

Note that the covering space for rays on a rf-dimensional torus is trivial 
since the momentum is single- valued. The hamiltonian in this case further- 
more depends on the modulus of the momentum in each direction and 
is positive definite, uq + V- = for a massless particle on the torus and 
cr n = for any periodic ray. The Casimir energy of a torus therefore is 
always negative, as was already found by deforming it and using Eq. (JHJ). 

4 Manifolds with Boundaries 

Estimates of the sign of the SCE are more difficult for manifolds with a 
boundary on which the scalar field satisfies some conditions. I here consider 
Dirichlet and Neumann boundary conditions only. Since classical rays reflect 
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specularly at a boundary, the basic strategy is to glue copies of the origi- 
nal manifold at the boundaries and consider the resulting covering manifold 
without boundary. At a boundary p = WS is discontinuous and there are 
(at least) two values for the momentum at any point near a boundary. By 

— » 

constructing the covering manifold for which p = VS is single-valued, one 
thus can treat a boundary in much the same way as a caustic. The phase 
retardation at a boundary depends on the boundary condition. The phase 
loss is 7r for Dirichlet and for Neumann boundary conditions. This ensures 
the correct behavior of semiclassical Green functions near the boundary. 

However, an additional correction to the semi-classical Casimir energy 
arises from periodic rays of the covering manifold that lie (entirely) within 
the boundary. The sign and magnitude of this additional contribution can 
be essential in determining the sign of the semi-classical Casimir energy and 
its dependence on the boundary condition. Let us examine some simple 
examples. 

4.1 Semiclassical Casimir Energy of a d-dimensional 
Half-Sphere 

Consider first the semi-classical Casimir energy of a half-sphere in (^-dimensions 
Space in this case is just cut in half at the equatorial Sj-i- Classical peri- 
odic rays in this case lie on two halves of great circles of the original S d with 
the same inclination that intersect on the equatorial Sj-i of the sphere. The 
momentum at any point on the half-sphere-annulus covered by a family of 
rays with fixed inclination therefore can take up to 4 values and one needs 

— * 

a 4-sheeted covering of this space to make VS single-valued. This cover- 
ing space is constructed in two steps. One first doubles the half-sphere and 
joins the two sheets at the equators to form an Sd- On this boundary-less 
double-covering of the half-sphere, geodesies are again great-circles and one 
again introduces two coverings for each annulus of Sd- Note that this last 
operation doubles the equatorial (boundary) Sj-i- One of these "equators" 
is where the upper and lower parts of the inner annulus join, the other is 
where the corresponding parts of the outer annulus join. This doubling of 
the boundary of the half-sphere cannot be avoided if the momentum on the 
covering space is to be single-valued. 

The periodic rays of this 4-sheeted covering of the ci-dimensional half- 
sphere evidently are those already found for the 2-sheeted covering space of 
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Sd- The only difference is that the phase of a ray may be retarded by n at 
every crossing of an "equator" (for Dirichlet boundary conditions). Since a 
periodic ray however, crosses the "equators" an even number of times, one is 
tempted to conclude that the semi-classical Casimir energy of the half-sphere 
is just half the semi-classical Casimir energy of Sd, irrespective of whether 
Neumann or Dirichlet boundary conditions have been imposed. 

However, this argument ignores classical periodic rays of the covering 
space that lie entirely within the boundary. The contribution of such rays 
in general depends on the imposed boundary condition. Since the field van- 
ishes on a boundary with Dirichlet's condition, classical periodic rays that 
lie entirely within the boundary should not contribute to Green's function 
and the spectral density. Their contribution to the spectral density of the 
manifold without boundary has to be subtracted in this case. 

Often, as in the case of a sphere, symmetry arguments can be invoked 
to relate the Casimir energy for Dirichlet and Neumann boundary condi- 
tions. Due to symmetry under reflections about the equatorial plane, eigen- 
functions of the hamiltonian can be chosen to satisfy either Neumann or 
Dirichlet boundary conditions at the equator of Sd- The sum of the Casimir 
energies on the half-sphere £(62/2; N) and £(62/2; D) for Neumann, respec- 
tively Dirichlet boundary conditions, therefore is the Casimir energy of the 
full sphere, 

£ c (S d /2; N) + £ c (S d /2; D) = £ c {S d ) . (19) 

The difference £ c (S d /2; N)—£ c (S d /2; D) is due to contributions to the Casimir 
energy from the boundary, i.e. due to periodic rays on the equator. The mag- 
nitude of this contribution in general is difficult to obtain without explicit 
calculation. It is not simply related to the Casimir energy of a d — 1 dimen- 
sional sphere, because £ c {Sd-i) does not include fluctuations transverse to 
the equator. One nevertheless can argue the sign of the difference, that is 
whether Neumann or Dirichlet boundary conditions lower the Casimir en- 
ergy. The point is that families of periodic rays in the vicinity of the equator 
of Sd are similar to those on Sd-i except that they nevertheless pass two 
caustics of order d — 1 - rather than of order d — 2 as for Sd-i - in every 
revolution. This gives an additional phase loss of n for every revolution no 
matter how close to the equatorial hyperplane the periodic rays are. Since 
the sign of the contribution of a periodic ray to the Casimir energy of S d -i 
is given by a(Sd-i) = (2n — l)(d — 2), the additional phase loss of nn on Sd 
changes the index to a(equator S d ) = (2n — l)(d — 2) + 2n. From this one 
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obtains that, 

{< 0, for d = mod 4 

= 0, for odd d . (20) 

> 0, for d = 2 mod 4 

The explicit calculation of the SCE of a 2-dimensional half-shell given in 
Eq. (J39j) of Appendix A is in agreement with the estimate of Eq. (J2(Jj) . It is 
well known that the Casimir energy is the same for a half-circle (i.e. interval) 
with Dirichlet, respectively Neumann boundary conditions at the endpoints 
(since there is no curvature for d — 1 this also coincides with ref. [10]). The 
Casimir energy due to a scalar on half an Einstein universe has also been cal- 
culated explicitly [34] and was found to be just half of the Casimir energy for 
the full Einstein universe irrespective of the boundary conditions (which cor- 
responds to d — 3 in Eq. (j2~Uj) . Note that the fate of the implicitly subtracted 
infinite vacuum energy proportional to the surface area of the equator of the 
halved Einstein universe was ignored in [33]. Fulling has pointed out;6. that 
(infinite) changes in the vacuum energy of the universe could be absorbed in 
the cosmological constant. They also may be cancelled by similar (infinite) 
contributions from other fields as in the case of super-symmetry. 

The following example of a (^-dimensional parallelepiped shows that the 
semiclassical evaluation of the Casimir energy of spaces with boundaries that 
are not smooth and intersect on lower dimensional manifolds can be even 
more involved. 

4.2 Semiclassical Casimir Energy of a d-dimensional 
Parallelepiped 

The Casimir energy of a massless scalar field confined to the interior of a 
parallelepiped with dimensions l\ x l 2 x . . . x 1$ has previously been obtained 
by Ambjorn and Wolfram |22j. They considered Neumann, Dirichlet as well 
as periodic and electromagnetic boundary conditions on the surface of the 
parallelepiped. I here give a more geometrical interpretation of some of their 
results in terms of periodic rays of a covering space. 

It suffices to consider the case where all the k of the parallelepiped are 
finite. The result of ref.[22j for a parallelepiped with some sides that are much 
longer than all others is found by taking the appropriate limits. Lacking a 
more concise notation, the SCE of a parallelepiped with dimensions l± x . . . x 
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will be denoted by, 

£ c (h, ■ ■ ■ , In] In+i, ■ ■ ■ i In+d] In+d+i, ■ ■ ■ , Id) ■ (21) 

Here Neumann boundary conditions are satisfied on < N < d pairs of 
parallel hypersurfaces that are distances h, . . . Jn apart, Dirichlet boundary 
conditions are satisfied onO<D<d — N pairs of parallel hyper-surfaces 
that are distances In+i, ■ ■ ■ Jn+d apart and periodic boundary conditions are 
assumed to hold on the remaining < d — N — D pairs of parallel surfaces. 

The symmetry of a parallelepiped implies that eigenfunctions satisfying 
periodic boundary conditions on a pair of faces, are even or odd under re- 
flection of the parallelepiped about these faces. This leads to the follow- 
ing relation between the Casimir energies of scalar fields satisfying different 
boundary conditions on the surfaces of d- dimensional parallelepipeds, 

£ c (h, ■ ■ ■ j In] In+i, ■ ■ ■ , In+d] %h, In+d+i, ■ ■ ■ > Id) = 

£ c (h, . . . , In] In+i, • • • , In+d] In+d+u ■ ■ ■ , Id) (22) 
+£c(h, ■ ■ ■ 5 In] h, In+i, ■ ■ ■ , In+d] In+d+i, ■ ■ ■ , Id) ■ 

Eq. (|2*2*|) is a relation between the Casimir energies of a scalar field on 
parallelepipeds where the periodic boundary conditions on a pair of parallel 
surfaces are replaced by Dirichlet or Neumann boundary conditions and the 
distance between the two surfaces is halfed. The total volume of the manifolds 
on the left- and right-hand sides of this equation thus are the same. The 
subtracted terms of the spectral densities are the same as well and the Casimir 
energies indeed are comparable. 

The spectrum of a parallelepiped differs only in a zero frequency mode 
for Neumann and Dirichlet conditions on a set of parallel surfaces. This fre- 
quency does not depend on the separation of the two surfaces. One therefore 
also has that, 

£ c (h, ■ ■ ■ 5 In] In+i, • • • j In+d] In+d+i, ■ ■ ■ , Id) — 

£ c (h, ■ ■ ■ , In] In+i, • • • j In+d] In+d+i, ■ ■ ■ , Id) (23) 
—£ c (h, ■ ■ ■ , In] h, In+i, ■ ■ ■ > In+d] In+d+i, ■ ■ ■ ,ld) ■ 

Note that the subtractions in the spectral density proportional to the d- 
dimensional volume cancel on the right hand side and the leading remaining 
subtraction is proportional to the volume of the (d — l)-dimensional paral- 
lelepiped on the left hand side. 
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Combining Eq. (|22|) and Eq. (|23|) one obtains the following recursive 
relation for the SCE of a parallelepiped, 



2£ c (li, . . . , l N ; Zjv+i) • • • 5 In+d] In+d+i, ■ ■ ■ ,h) — 

£c{h, ■ ■ ■ , In', In+i, ■ ■ ■ , In+d', 2Zi , In+d+i, ■ ■ ■ Jd) 
+£ c (h, ■ ■ ■ ,In] In+i, ■ ■ ■ , In+d', In+d+i, ■ ■ ■ , Id) 



(24) 




, In; Zi, In+i, ■ ■ ■ 5 In+d; In+d+i, ■ ■ ■ Jd) = 
-c{h, ■ ■ ■ ,In] In+i, ■ ■ ■ j In+d] %h, In+d+i, ■ ■ ■ ,ld) 
-£c(h, ■ ■ ■ ,In] In+i, ■ ■ ■ , In+d', In+d+i, ■ ■ ■ , Id) ■ 



Eq. (J24)) expresses the Casimir energy of a parallelepiped in terms of 
Casimir energies of parallelepipeds with non-periodic boundary conditions 
on fewer sets of parallel plates. Repeated application of these relations thus 
gives the Casimir energy of a parallelepiped with Neumann, Dirichlet and 
periodic boundary conditions in terms of the Casimir energies of tori only. 
Thus the Casimir energy of a three-dimensional parallelepiped with Dirichlet 
boundary conditions on all six faces may be decomposed as, 



The corrections proportional to the Casimir energies of lower dimensional 
tori are due to periodic rays of the covering space that lie on the boundaries 
of the original parallelepiped. Since the subtraction for a torus is the Weyl- 
contribution proportional to its surface, the subtractions for a general paral- 
lelepiped include Weyl-terms proportional to the "area" of its hyper-surfaces, 
"lengths" of the intersections of its hyper- surfaces etc.. The Casimir energy 
of a general parallelepiped thus can be compared with that of other systems 
of the same volume, area of the boundary (with the same boundary condi- 
tions), length of intersections of hyper-surfaces etc... The simplest class of 
spaces that satisfy all these conditions are c?-dimensional generalizations of 



S£ c (; li, 1%, l^; ) 



4£ c (;l 2 ,l 3 ;2h) - A£ c (;l 2 ,l 3 ;) 
2£ c (;h;2l 1 ,2l 2 )-2£ c (;l 3 ;2l 1 ) 
-2£ c (;l 3 ;2l 2 )+2£ c (;l 3 ;) 
£ c (;;2h,2l 2 ,2l 3 ) 

—£ c (; ; 2li, 2l 2 ) — £ c {; ; 2/ 1; 2/ 3 ) — £ c {; ; 2l 2 , 2l 3 ) 
+£ c (;;2h)+£ c (;;2l 2 )+£ c (;;2l 3 ) . 



(25) 
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Power's box^J] with a fixed number of orthogonal but movable walls 4 

Let us turn to the construction of the appropriate covering space for a 
parallelepiped. The momentum is single valued only when periodic boundary 
conditions hold on all pairs of hyper- surfaces. N + D = and there is no 
need to introduce additional sheets. The result is the same as for the d- 
dimensional torus in section 12.11 the Casimir energy in this case is negative 
for any dimension of the torus and any lengths, Zj, of its cycles. 

For N + D > 0, the momentum is not single valued. Each pair of faces 
with non-periodic boundary conditions requires a double covering since the 
component of momentum that is perpendicular to these surfaces can have 
either sign. One recursively constructs this covering space as follows. 

Consider first the pair of faces with coordinate x\ = and x\ — l\. The 
boundary condition on this pair is not periodic and the first component of 
momentum (for rays of fixed energy) therefore is double valued. It is single 
valued on a covering space obtained by joining a second sheet of the original 
parallelepiped to the first at the boundaries X\ — and x\ = l\ to form 
a cylinder-like covering space. Periodic rays that reflect from the X\ — 
and x\ = l\ faces of the original parallelepiped pass smoothly through these 
borders from one sheet to the other in the covering space. Although one must 
keep track of the position of the original boundaries at x\ = and x\ = l\, 
the problem of constructing a covering space on which momentum is unique 
has been reduced to that of a parallelepiped with only D + N — 1 pairs of 
hyper-surfaces on which non-periodic boundary conditions hold. Note that 
the first dimension of this covering space is now twice that of the original 
parallelepiped. 

li N + D > 1, the procedure is repeated with the pair of faces at X2 = 
and 22 — h of this double cover of the original parallelepiped. This results 
in a covering parallelepiped with non-periodic conditions in just iV + D — 2 
dimensions. The length of this covering parallelepiped in the second dimen- 
sion is again twice that of the original parallelepiped. The process ends with 
the pair of faces at xd+n = and xd+n = Id+n- 

The covering space of a general d- dimensional parallelepiped thus is a 
d-dimensional torus with cycles 2Z 1; . . . , 2l D+N , l D+N+1 , ...,1a- The original 

4 Power's original construction of a box with just one movable wall for the Casimir force 
between two parallel plates has been reexamined and extended in|32|. Since the surface 
areas of the two parallelepipeds are linearly dependent on their volume, considering a box 
with just one movable wall does not isolate the surface dependence of the Casimir energy 
of a parallelepiped - one has to consider more than one movable wall. 
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parallelepiped is covered 2 D+N times. The 2(D + N) (d — 1) -dimensional 
hyper-surfaces with non-periodic boundary conditions of the original paral- 
lelepiped now are (d — l)-dimensional interfaces between sheets of this cov- 
ering space at x\ — 0, X\ = h,X2 = 0, . . . , xd+n = and xd+n = Id+n- 

Since a general periodic ray crosses pairs of boundary surfaces, /3 7 = 
for periodic- and Neumann, respectively Dirichlet boundary conditions on 
opposing pairs of hyper-surfaces 5 . 

On this toroidal covering space there are periodic rays that lie entirely 
within the d — 1-dimensional hyper-surfaces that are projected onto the 
boundaries of the original parallelepiped. The contribution to the Casimir 
energy due to these rays depends on the imposed boundary conditions. The 
hyper-surfaces of a parallelepiped in addition intersect on lower dimensional 
hyper-edges that also contain periodic rays of the covering space. 

The corrections due to rays on these lower-dimensional tori are clearly 
visible in Eq. (J23jl and Eq. (|25|). They are in one-to-one correspondence with 
the boundary surfaces, edges, etc. of the parallelepiped. The contribution 
of rays on the lower- dimensional surfaces has to be subtracted for Dirichlet 
boundary conditions because the field vanishes on this surface of the paral- 
lelepiped in this case. That the lower-dimensional correction in Eq. (J2%jl is 
added for Neumann conditions then follows from the reflection symmetry of 
the parallelepiped, which implies Eq. (J22J). 

The sign of the Casimir energy of a parallelepiped in general therefore 
is determined by the boundary conditions on its surfaces. One can argue 
that the Casimir energy due to periodic rays of a boundary surface is always 
negative when this surface is embedded in a higher-dimensional space of 
vanishing curvature: there are no caustics to contend with and the energy 
surface has no zero modes. 

Whether the Casimir energy due to periodic rays on such a boundary 
surface has to be subtracted or added depends on the imposed condition. The 
sign of the overall Casimir energy of a parallelepiped therefore depends on the 
relative magnitude of boundary contributions with opposite sign. Which sign 
prevails in general will depend on the actual dimensions of the parallelepiped. 
The previous analysis nevertheless allows for a few general statements about 

5 /3 7 = even > 2 if the scalar field satisfies Neumann's condition on one, but Dirichlct's 
condition on the other of a pair of parallel hyper-surfaces. Such asymmetric boundary 
conditions on a pair of parallel surfaces can give a change in sign of the semi-classical 
Casimir energy [Tlll35j . It is due to a phase loss of tt for some dominant primitive periodic 
rays and can be explained along the lines used for manifolds without boundarv|14|. 
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the sign of the Casimir energy of a parallelepiped: 

• Since contributions from periodic rays on lower dimensional boundary 
surfaces are always negative, the sign of the SCE of a parallelepiped 
with only Neumann and periodic boundary conditions is negative in 
any dimension. 

• Replacing Dirichlet- by Neumann- boundary conditions on a pair of 
parallel surfaces decreases the Casimir energy of the parallelepiped. 

• If one dimension of the parallelepiped is much smaller than all others, as 
for two parallel infinite hyper-planes, the difference in Casimir energy 
for Neumann and Dirichlet boundary conditions tends to vanish (since 
all contributions to the Casimir energy due to rays on lower dimensional 
surfaces become negligible). 

5 Conclusion 

By defining the SCE through the required subtraction po(E) in the spectral 
density of a system one describes a class of systems whose vacuum energies 
can be compared. By definition the subtracted spectral density p(E) of 
Eq. (jUJ) gives a finite Casimir energy. It is approximated semiclassically by 
contributions due to classical periodic rays. For a massless scalar field on d- 
dimensional spheres and tori as well as on related spaces with boundary such 
as half-spheres and parallelepipeds, this SCE coincides with other definitions 
whenever the subtractions in the spectral density are the same for a non- 
trivial class of systems(see Appendix A). 

[One may argue that certain systems, such wedges formed by two semi- 
infinite planes joined at the common edge [3J, do not have classical periodic 
rays and that this semiclassical approach may thus be rather limited in scope. 
However, the implicit subtraction is not universal in this case, and depends 
on the opening angle of the wedge in a non-linear fashion that prohibits 
comparisons between systems with wedges of different opening angle. The 
extracted finite Casimir energy of a particular wedge in this case is of little 
physical significance and for instance does not determine the torque between 
the two plates of the wedge.] 

The geometrical description of the SCE in terms of periodic rays gives 
insights into qualitative features of this part of the vacuum energy that oth- 
erwise are rather mysterious. The need for an explanation of the sign of 
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Casimir energies was very nicely formulated by J. Sucher: "Understanding 
the signs is a sign of understanding" [35] . I have here presented some evi- 
dence that the sign of the SCE depends critically on optical properties of 
"important" (short) periodic rays. 

The semiclassical contribution to the Casimir energy due to an isolated 
periodic ray has a definite sign (see Eq. ©). The contribution due to a 
class of periodic rays of an integrable system also is of a definite sign (see 
Eq. ()16|). For isolated periodic rays, the sign is determined by the winding 
number cx 7 of their stable and unstable manifolds .25 , whereas it essentially 
is given by the Keller-Maslov index f3 n of a class of periodic rays in integrable 
systems. [When the Hessian of Eq. (114)) is not positive definite, the sign of 
the contribution of a class of rays is more appropriately given by Eq. (|T7|) .] 

The semiclassical expressions of Eq. © and Eq. (|TB|) suggest that the 
Casimir energy very often is dominated by the shortest periodic rays or class 
of periodic rays. Arguing that the contribution of the shortest periodic rays 
is continuous under small deformations of the manifold one can estimate 
the sign of the Casimir energy of integrable systems in two ways: either by 
direct computation of the index in Eq. (|17|) . or by slight deformation of the 
integrable system and computation of the index of Eq. ((7|) for the shortest 
periodic rays. Continuity under deformations also explains the pattern of 
signs for the Casimir energies of a massless scalar on spheres and of tori 
of various dimensions as well as changes in sign when spheres are strongly 
deformed. Sign changes of the SCE in these cases are accompanied by a 
change the number of conjugate points of the shortest rays. 

The sign of the SCE is harder to determine for spaces with boundary. 
The shortest periodic rays again dominate, but the sign (and magnitude) of 
contributions due to periodic rays that lie within the boundary depends on 
the boundary conditions. Such boundary rays can be among the shortest 
periodic rays and in some cases dominate the Casimir energy. 

Only (^-dimensional half-spheres and general parallelepipeds with Neu- 
mann and Dirichlet boundary conditions were considered. However, the gen- 
eral arguments remain valid in more realistic situations with, for instance, 
electromagnetic fields. The sign of the Casimir energy of spherical- and 
cylindrical- cavities with idealized metallic boundary conditions can appar- 
ently be understood in terms of the phases of the shortest periodic ravs|36j. 

Determining the sign of a SCE in this sense is reduced to a problem of 
geometrical optics. For certain simple manifolds, such as tori and spheres, one 
finds that all contributions due to periodic rays are of the same sign. The sign 
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of the SCE is unambiguous in this case. The fact that classical dynamics of 
periodic trajectories seems to determine the sign of the Casimir energy raises 
the intriguing question whether the sign is a topological characteristic of the 
phase space. In general, and in particular for manifolds with boundaries, 
periodic rays of comparable length contribute to the SCE with opposite sign. 
It then depends on the boundary conditions and metric characteristics of the 
space whether the SCE is positive or negative. 
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A Casimir Energies and Curvature Correc- 
tions: Spheres and Half-Shells 



I here explicitly compute the Casimir energies due to a scalar field for low- 
dimensional spheres Si, S 2 , S 3 and S4 as well as for the two-dimensional half- 
shell with Neumann, respectively Dirichlet boundary conditions. Of special 
interest are the associated subtractions in the spectral density that render 
the Casimir energies finite. As explained in the main text, these subtractions 
determine classes of systems with finite vacuum energy differences. The 
finite Casimir energy is given (exactly) by contributions due to periodic rays 
whenever the subtractions are universal. 

A.l The circle 

The Casimir energy of a scalar field on Si probably is the most transparent 
example. The curvature vanishes and the energy spectrum of a massless field 
on Si clearly is Ei = IHc/R for integer / > 0. For I > the energy eigenvalues 
are 2-fold degenerate. The vacuum energy of a massless scalar on Si thus is 
formally given by, 




(26) 



1=1 
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One may regularize this divergent expression in many ways. One of the 
more popular is by analytic continuation in the exponent s of the for s < 
— 1 manifestly convergent sum ((— s) = Y1T^ S - This method, known as 
zeta function regularization, has been claimed[30 to be not just the most 
elegant, but also the only rigorous mathematical definition of the sum in 
Eq. (J26j) . That may be true but gives appreciably little insight into the 
physical significance of the finite Casimir energy one obtains, which for a 
circle is — j^. 

Let us subtract from the spectral density ps ± (E; Si) = J^'ilo 25(E—hcl/R) 
the smooth Weyl density po(E) = The SCE of a circle then is defined 
to be, 
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(27) 



The second line of Eq. (|2*7]l expresses the subtracted spectral density in terms 
of the semiclassical contribution due to periodic rays of length 2nRn, n > 0. 
The final answer coincides with that of zeta function regularization, but one 
now explicitly knows the implicit subtraction in the spectral density required 
to obtain it. The subtracted Weyl contribution to the spectral density is 
proportional to the circumference of the circle. The difference in Casimir 
energies of two circles thus is not the difference in their vacuum energy. It 
in fact would cost an infinite amount of energy to change the radius of the 
circle by a finite amount. 

The subtraction nevertheless is universal in the sense that the (so defined) 
finite Casimir energy reproduces all derivatives of the vacuum energy with 
respect to the radius of 2 nd and higher order. The difference in total Casimir 
energy for instance gives the energy required to change the relative radii of 
disjoint circles while keeping their sum constant 6 . 

6 The situation is unstable due to the minus sign of the Casimir energy in Eq. (|27|l . It 
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A.2 The Two-Sphere 

The 2-dimensional sphere S 2 is a prototypical manifold with positive curva- 
ture. The Casimir energy of a scalar field on S 2 depends on how the scalar 
couples to the curvature. It is known [T%1 |2~9"1 Id] . and we shall see, that this 
coupling to the curvature is of a particular strength for a massless field. 
However, to better compare with ref. [3U], let us for the moment ignore the 
coupling to the curvature. 

The eigenf unctions of the Laplace-Beltrami operator on a two-sphere are 
the spherical harmonics with (21 + l)-fold degenerate eigenvalues. Without 
coupling to the gaussian curvature of S 2 , the energy eigenvalues of a scalar 
field are Ei = Hcy/l(l + l)/R. The vacuum energy of a two-sphere of radius 
R in this case is formally given by the divergent sum 

k 00 

£ V ac(S 2 ) = 2! E( 2/ + !) V^+l). (28) 

1=0 

Zeta function regularization (ZR) of Eq. (J2~%|) gives [30^ the negative value, 

£zr(S 2 ) = Km ^ f> + + = -0.132548^ (29) 
s ^ 1=0 

for this Casimir energy. A physically perhaps more transparent treatment of 
the sum in Eq. (}28|) by heat kernel regularization (HK) gives, 

£ HK {S 2 ) = hm ^ jf> + l)VKlTT)e-^ - ^ J = -0.13254 8 |. 

(30) 

Subtracting the integral over I £ [0, 00} from the sum over angular momenta 
thus leads to the same finite result for the Casimir energy as zeta function 
regularization. However, no obvious physical principle apart from finiteness 
of the result seems to dictate the particular subtraction in Eq. (}30|) . [For a 
justification and the relation to subtractions in other regularization schemes 
see ref. 

A straightforward physical interpretation of the subtraction is possible in 
the semiclassical treatment: one again subtracts a smooth "classical" part of 
the spectral density. po(E) should not depend on the detailed shape of the 

implies that the energy is minimized by shrinking all but one of the circles to a point. 
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surface. For a 2- dimensional manifold without boundaries, one may subtract 
the classical Weyl contribution to the spectral density proportional to the 
area A of the manifold, po{E) = AE/(2Tr(hc) 2 ). The subtracted spectral 
density is, 



p{E)=p{E)-p {E) 



oo 
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p(-E') can be expressed in terms of semiclassical contributions from periodic 
rays PI EZ| , 
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One verifies that Eq. (|32Jl is equivalent to Eq. (}3T]) by recognizing that apart 
from an = fl term, the real part of the sum over n in Eq. (}32j) is a periodic 
^-distribution. The sum in Eq. (f3*2*|) on the other hand can be interpreted as 
due to contributions from classical periodic rays that wind n times about a 
geodesic of the two-sphere with a total length L n = 2nRn. 

The Casimir energy corresponding to the spectral density p(E) of Eq. (}3*2*|) 



is. 
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where the last number was obtained by numerical evaluation of the integral 7 . 
The SCE thus coincides with that of zeta function- and heat kernel- regular- 
ization. The Casimir energy of a spherical shell is completely described in 



Substitution of the variable i£(E) = y 1 + (^^r) — 1 an d a clockwise rotation of the 
integration contour by 90° gives the integral that converges exponentially and is numeri- 
cally well behaved. 
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terms of its classical periodic rays. Moreover, this procedure gives a trans- 
parent physical meaning to the subtraction: the spectral density p(E) is the 
difference to the universal Weyl term in the spectral density for 2-dimensional 
surfaces of the same area. The Casimir energy thus does not give the en- 
ergy required to change the radius of S2' the subtracted (divergent) term 
proportional to the surface area of the sphere is far more important for this 
energy difference (22] ■ The Casimir energy on the other hand does allow 
the computation of the energy required to deform the sphere into another 
smooth manifold without boundary of the same total area. It thus for in- 
stance makes sense to speak of the energy required to change a two-sphere 
to a discus or an elongated ellipsoid of the same area. The energy required 
for such deformations is finite. 

The effective action of a periodic ray of arc length L n = 2irRn in Eq. (}3*2*j) 

is, 

S n = jv&* = p(E)L n} (34) 

with 

p 2 {E) = {E/c) 2 + h 2 /{AR 2 ) . (35) 

Eq. ()35j) is not the dispersion relation one expects for a massless particle. 
p 2 {E) differs from (E/c) 2 by a term of order h 2 that is proportional to the 
Gaussian curvature k = 1/R 2 of S2. The curvature of the manifold results 
in a dispersion relation that would corresponds to a tachyon with velocity 
v(E) = (dp/dE)' 1 = cy/l + (hc/(2RE)) 2 > c. This can also be seen by 
rewriting, 

Interpreting H(l + 1/2) jR as the eigenvalue of the momentum operator for 
a field satisfying anti-periodic boundary conditions on a circle of radius R, 
this dispersion clearly is tachyonic: the effective mass squared is negative 
m 2 = —(h/2cR) 2 . The fact that the effective action in Eq. ()34j) does not 
vanish for E = is another indication that the spectrum of the Laplace- 
Beltrami operator on curved spaces is not the spectrum of a massless scalar 
field. 

A particle is massless if its dispersion is p(E) = E/c. The generalization of 
the wave equation for <i-dimensional Riemannian manifolds A4 with metric 
gij and Gaussian curvature k includes a coupling to the curvature. The 
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appropriate wave equation for a d + 1-dimensional (ultrastatic) space-time 
with curvature is = c~ 2 d 2 <p, where 

d ■■ d (d-1) 2 
A.HVVa^^-^- (37) 

This curvature correction to the Laplace-Beltrami operator of M. also arises 
naturally from the measure of the path integral and is consistent with the 
required short-time behavior of the Feynman propagator [21]. The particular 
strength of the coupling to the curvature in Eq. (jH7|1 preserves conformal 
invariance 8 of the wave equation[6j. For a <i-dimensional sphere Sd of radius 
R, k = 1/R 2 and the eigenvalues of — A d defined in Eq. (J3~7|) are (I + (d — 
l)/2) 2 /R 2 for integer I > 0. The Casimir energy due to a massless scalar 
field on a two-sphere thus is, 



S C (S 2 ; m = 0) 
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1 AtcR 2 ^ r°° r / 2ttR \ 

= Im r75E^ z r : = - ( 38 ) 

4tt 3 R ^ n 3 

n=l 

That the Casimir energy of a massless scalar field vanishes for a two-sphere 
agrees with the arguments presented in Sections 12.11 and El Zeta function- 
and heat kernel-regularization also give a vanishing Casimir energy when the 
coupling of a massless scalar to the curvature is included. Note that the 
average spectral density of a 2-dimensional manifold of the same area again 
was subtracted to obtain the finite answer of Eq. (jSHj)- Although the second 
line of Eq. ()38|) is a sum over contributions due to periodic rays of the two- 
sphere, it does not have the form of Eq. The amplitude in particular is 
proportional to the energy here. As discussed in Section El free motion on a 
two-sphere is integrable and the corresponding periodic rays are not isolated. 
As shown in Appendix B, a semiclassical treatment along the lines of ref. j 



8 k in Eq. I|37|l is the Gaussian curvature of M. rather than the Ricci curvature scalar 1Z 
of the space-time RxM considered injSJ. For an ultrastatic space-time the two curvatures 
are related by 1Z — k(<1— l)/d and the particular coupling strength in Eq. Ij37|l corresponds 
to the conformal coupling £ = + — 2)/(4((g?+ 1) — 1)) discussed injS]. 
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gives the semiclassical response function of Eq. (|38|) without reference to the 
exact quantum mechanical spectrum. 

Let us now turn to a manifold with boundary and consider a two-dimensional 
half-sphere S 2 /2 with Neumann (N), respectively Dirichlet (D) boundary 
conditions on the equatorial circle 9 Due to the symmetry of the two-sphere 
upon reflection about its equator, eigenf unctions of the Laplace-Beltrami op- 
erator either satisfy Neumann or Dirichlet conditions on the equator and 
thus £(S2/2;N) + £(S2/2;D) = £(82) = 0. Since one more mode satisfies 
Neumann's boundary condition than Dirichlet's for every energy eigenvalue, 
one has for the Casimir energy of a half-sphere, 
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This again agrees with zeta function regularization for this case. Note that 
an additional universal subtraction is necessary. It is a Weyl contribution to 
the spectral density proportional to the length of the boundary. Due to the 
subtraction of terms proportional to the area of the surface and proportional 
to the length of the boundary, the Casimir energy of a massless scalar on 
the half-sphere can be compared to the Casimir energy on another smooth 
2-dimensional manifold of the same area and with a smooth boundary of the 
same length only. 

Note that the necessary subtraction for a half-sphere with boundary is 
not universal if the curvature correction to the Laplace-Beltrami operator 
in Eq. ()37|) is ignored. In this case one obtains the following difference in 
Casimir energies, 



£(S 2 /2;N) - £(S 2 /2-D) = 

9 A more general treatment of the Selberg trace formula for two-dimensional manifolds 
with boundaries is given in^Uj. However, the particular example of a massless scalar on 
a half-sphere, does not appear to satisfy the restriction that h(p) = e^ p ^ be analytic on 
the strip |Imp| < 1/2 + e. 
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The subtracted spectral density po(E) in Eq. PU|) is again proportional to the 
length of the boundary, but it now depends on the Gaussian curvature k = 
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1/i? 2 as well. Expanding for small values of the curvature Ej y E 2 + 

l-/t(nc) 2 /(8 J B 2 ) shows the presence of a (logarithmic) divergent contribution 
to the vacuum energy that is proportional to 1/R. No derivative of the 
vacuum energy with respect to the radius of the half-shell is finite and the 
subtracted (finite) Casimir energy of Eq. (jlUj) is of no physical significance. 

In zeta function regularization the logarithmic divergence manifests itself 
as a pole. The principal value prescription ignores this pole contribution 
to the Casimir energy to produce a finite value[30j of — 0.166080/lc/i? for 
the above difference of the Casimir energies of half-shells with Neumann 
and Dirichlet boundary conditions. Since the implicit subtraction is not 
universal, it should not surprise that the value ascribed by principal value 
prescription differs from the one in Eq. ()4Uj) - obtained by prescribing a 
particular subtraction of the spectral density. The finite part of the Casimir 
energy of this system is quite arbitrary since its vacuum energy cannot be 
compared to that of any other system anyway. The appearance of a pole in 
zeta function regularization in this sense indicates that the associated implicit 
subtraction cannot be universal and that the finite value obtained by some 
prescription is not very meaningful. 

If one is of the opinion that a Casimir energy due to a physical field should 
enjoy some degree of universality, small deformations of a smooth manifold 
should not require an infinite amount of energy. The lack of universality of 
the Casimir energy on a half-shell in this sense is another indication that the 
spectrum of the Laplace-Beltrami operator is not that of a physical particle. 
This favors Eq. (|37|) as the appropriate generalization of the wave operator 
to manifolds with curvature. The argument perhaps is more convincing if 
one notes that even the Casimir energy of a scalar on S3 (a manifold without 
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boundary that could serve as a model of the spatial part of our universe) is 
not universal [SU] without curvature correction. 



A.3 S3 and S 4 

For comparison with the results of ref.j3U] and in support of the arguments 
for the sign of the Casimir energy of a massless scalar field on d-dimensional 
spheres of Sections 12.11 and El I here give the universal Casimir energies of a 
massless scalar field on d = 3 and d = 4-dimensional spheres. I have included 
the appropriate coupling to the curvature of Eq. (}3Tj) and the results therefore 
differ from those without curvature correction injHU]. The difference is not 
minor: the Casimir energy of a scalar on S3 in particular is not universal and 
devoid of physical implications without coupling to the curvature. In zeta 
function regularization this is indicated by the appearance of a pole|3Uj. 

The Casimir energies due to a scalar on S3 and S4 can be obtained in 
much the same manner as for the two-sphere - the main difference is the 
degeneracy of the spectrum of the Laplace-Beltrami operator in Eq. (|37|) . It 
is (I + l) 2 for S 3 and (Z + 1)(Z + 2) (21 + 3)/6 for S 4 . As can be seen from 
Eq. ([37)1 . the spectrum of a massless field with the appropriate coupling to 
the curvature is E\ = Hc(l + 1)/R on S3 and E\ = hc(l + 3/2) /R on S4. One 
then finds, 

Tie fxc 
S c (S 3] m = 0) = — = +0.0041666 ... - 

S C (S 4 ; m = 0) = . (41) 

These results again are in agreement with those of zeta function regulariza- 
tion for this case. There are no logarithmic divergent subtractions and no 
poles appear in zeta function regularization when the appropriate curvature 
correction is taken into account. For S 4 the subtraction of the Weyl contri- 
bution proportional to the 4-volume of S4 is not sufficient to obtain a finite 
Casimir energy. An additional smooth spectral density proportional to R 2 
must be subtracted and, 
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The additional term in the smooth part of the spectral density is propor- 
tional to the integrated curvature of the manifold. One thus can compare 
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the (vanishing) Casimir energy of a massless scalar field on S4 with that of 
other 4-dimensional manifolds without boundary of the same J^-volume and 
average curvature. 

B Semiclassical Derivation of the Spectral Den- 
sity of a Scalar on a Two-Sphere 

For completeness I here give the semiclassical calculation of the Casimir en- 
ergy ^ C (S , 2 ) of Eq. (|3*H|) for a massless scalar particle on a two- sphere without 
reference to the quantum mechanical spectrum. Although a somewhat triv- 
ial example of the more general asymptotic expansion for elliptic self-adjoint 
pseudo-differential operators on Riemannian manifolds [3H1 EH]- it does illus- 
trate some of the basic ideas at an elementary level and illustrates the general 
procedure of section El 

The system has two constants of motion in involution and is thus inte- 
grable. One of these can be taken to be the square or magnitude of the 
momentum \p\, the other the angular momentum l z . The corresponding ac- 
tion variables for a two-sphere of radius R are ii = R\p\ > and I2 = l z with 
|-?2 1 < I\- The hamiltonian for a free massless particle on the two-sphere is 
iy(Ji,/ 2 ) — c \p\ — cl\jR. Following ref . (2H1 EE] > the oscillating part of the 
semiclassical spectral density is, 

1 00 ' poo plx 

p{E; S 2 ) = - / / dI ^ E - F ( J i> J 2)) e 2 ^ nh+mh)/h - n/2] 

n n,m=-oo Jo 

(43) 

The integration domain of the action variables in Eq. (}4*3*|) is due to classical 
considerations only and the primed sum in Eq. (J43j) here implies that the 
(classical) contribution with m = n = is to be omitted. Note that the 
m = n = 0-term is just the po{E) that is subtracted from the full spectral 
density. The spectral density evidently is given in terms of contributions 
from periodic trajectories that wind (n, m) times about the (in this case two- 
dimensional) invariant tori of constant energy. The phase retardation by nir 
of the family of trajectories that wind n-times about the whole two-sphere is 
related to the uniqueness of the quantum mechanical wave function|31j and 
might not be expected classically. Keller's construction in section l3~T1 shows 
that this phase loss is due to the two caustics of first order formed by families 
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of periodic rays of the same inclination 10 . Inserting Eq. (J43j) in Eq. (@J) and 
performing the integral over E one obtains, 



i 00 1 poo rli 

£c(S 2 ) = — Yl J dhj ^dI 2 H(h,I 2 )e 2m ^ h+mI ^ h - n W . (44) 



Although a direct consequence of Eq. (}4*3*j) . it perhaps is remarkable that the 
SCE of an integrable system is expressible in terms of the Fourier-transform 
of the classical hamiltonian with respect to action variables. Since the hamil- 
tonian H(I\, I 2 ) = cl\/ R does not depend on I 2 and one is interested in the 
term of order h only, one can perform the integral over I 2 in saddle point 
approximation. This gives, 

f 1 dl 2 e 2nimh / h = h [ ^ dx e 2mmx = 2h8 m0 + 0(H) . (45) 

J-Il J 

The only sizable contribution to the I 2 integral thus is from the m = term 
of the sum. Inserting Eq. (f43|) in Eq. (jUj) one finally obtains for the SCE 
due to a massless scalar on 5*2, 

4?r 3 )R^ n 3 K ' 

The integral over Ii in Eq. (|46j) coincides with that over the energy in Eq. (}3*Hj) 
if one recalls that Ii = ER/c. Note that higher moments of the semiclassi- 
cal spectral density do not vanish. One in particular finds that the second 
moment of the semiclassical spectral density is negative, 

The contribution of periodic rays to the spectral density of a massless scalar 
thus cannot be positive semi-definite for S 2 . Upon subtracting the "classical" 
Weyl contribution to the spectral density there in fact is no reason why the 
remainder must be positive. The negative sign in Eq. (}4Tj) can be traced to 
the phase loss of ir for every revolution of a periodic ray and thus ultimately 
to the Keller-Maslov index of a class of periodic rays. 

10 For fixed energy the inclination of a periodic orbit is determined by l z = 1% on S2 
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